Melting artificial spin ice 
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Artificial spin ice arrays of micromagnetic islands are a means of engineering additional energy 
scales and frustration into magnetic materials. Here we demonstrate a magnetic phase transition in 
an artificial square spin ice and use the symmetry of the lattice to verify the presence of excitations 
far below the ordering temperature. We do this by measuring the temperature dependent mag- 
netization in different principal directions and comparing with simulations of idealized statistical 
mechanical models. Our results confirm a dynamical pre- melting of the artificial spin ice structure at 
a temperature well below the intrinsic ordering temperature of the island material. We thus create 
a spin ice array that has real thermal dynamics of the artificial spins over an extended temperature 
range. 
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Geometric frustration is observed in many physical sys- 
tems. The textbook example is the frustration of proton 
interactions in water ice, giving rise to proton disorder, 
as revealed in the pioneering experimental work of Gi- 
auque and Stout [ ] and the theoretical interpretation of 
Pauling [2] . Frustration in antiferromagnets analogous to 
the ice model was predicted theoretically by Anderson in 
1956 [3]. Only much later was ice- type disorder observed 
in magnets, and most surprisingly, this was in ferromag- 
netic materials that were thus named 'spin ice' [ - ]. The 
spin ice phenomenon relies on dipole-dipole interactions 
and the concept was later generalized to include arrays 
of magnetic islands with a spin ice type geometry, the ar- 
tificial spin ices [7-13]. Beyond the ice- type systems, the 
property of frustration is responsible for the occurrence 
of thermodynamically metastable phases in a variety of 
systems, including structural (amorphous materials) [14], 
magnetic (spin glasses) [1^^] and polymeric systems [16]. 
In general, the absence of unique ground states [ ] and 
the presence of kinetic constraints [17] makes the study of 
excitations extremely important. These, in combination 
with the degree of frustration, are the driving forces for 
disorder in such systems [ib]. An order-disorder transi- 
tion has been proposed for colloid spin ice model systems 
with adjustable interactions [19] and dynamical studies 
have already proven the importance of distortions [20]. 

The aim of the present work is to realize a thermally in- 
duced transition in an artificial spin ice array. In previous 
work on artificial magnetic spin ice, dynamics have been 
induced by changes in an applied magnetic field, rather 
than by thermal excitations [ - ]. The use of field- 
cycling protocols has enabled the realization of quasi- 
degenerate spin states reminiscent of real spin ice [ , ] , 
as well as ordered states that are in static equilibrium 
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FIG. 1. (a) Atomic Force Microscopy (AFM) image of the 
artificial square spin ice array. The major symmetry axes are 
indicated in the lower left corner, (b) Remanent magnetic 
configuration of the spin ice array after applying a saturation 
field parallel to the [10] direction and (c) parallel to the [11] 
direction, (d) The 16 possible remanent magnetic configura- 
tions for the artificial square spin ice vertices. After applying 
a saturation field along the [10] direction there is a fourfold de- 
generacy of the possible remanent vertex configurations. In 
comparison, after removing the [11] field, there is only one 
possible remanent vertex configuration. 



[ ]. All these results are obtained at temperatures far 
below the ordering temperature of the patterned struc- 
tures. Recent work by Morgan et at. [13] explored for 
the first time ground states in artificial spin ice, obtained 
during the growth of the pattern. 

Here we used J-doped Pd(Fe) to enable the exploration 
of the thermally driven order-disorder transition of the 
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FIG. 2. Normalized remanent magnetization versus temper- 
ature curves after a magnetic field has been applied parallel 
to the [10] and [11] directions, compared with the remanent 
magnetization of the continuous film used for patterning the 
arrays. The collapse of the array magnetization at a temper- 
ature well below the Curie temperature of the material (Tc 
= 230 K) is consistent with the appearance of thermally in- 
duced dynamics of the macro-spins that comprise the array. 
The inset shows representative normalized magnetic hystere- 
sis loops for a magnetic field applied parallel to the [10] and 
[11] directions, at a temperature of 12 K. 



artificial spin ice structure. The Curie temperature of 
(5-doped Pd(Fe) films can be adjusted by variation of the 
iron (Fe) layer thickness in the palladium (Pd) [25, 2i ]. 
A 1.2 monolayer thick Fe layer was embedded between 10 
monolayers of Pd, which resulted in a Curie temperature 
of Tc = 230 K. The film was patterned by an electron 
beam lithography and Ar ion milling into an artificial 
square spin ice array, of 750x250 nm^ islands with a 
period of 1000 nm, extending over a 1.5x1.5 mm^ area 
(Fig. 1). The total magnetic moment for one island was 
determined to be Mq = 1.1 x 10"^^ Am^ at 5 K. 

Magnetic hysteresis loops were recorded using the 
magneto-optical Kerr effect (MOKE) at temperatures 
from 5 K to 300 K [ ] (see inset in Fig. 2). The field 
was applied in two directions: (i) parallel to the [10] and 
(ii) parallel to the [11] directions of the patterned array, 
see Fig. 1. When the field is removed, different popu- 
lations of magnetic states for the different directions are 
obtained. For example, in the absence of thermal excita- 
tions, such a field cycling effectively removes vertices of 
type I, due to symmetry reasons. The normalized rema- 
nent magnetization, Mrem, of the magnetic island array 
is compared to the continuous sheet material in Fig. 2. 
The remanent magnetizations follow a similar trend up 
to 100 K, but at higher temperatures the magnetic is- 
land array crosses over to a regime in which the magne- 
tization is strongly reduced, approaching zero at around 
200 K, while that of the continuous film falls to zero at 



Tc = 230 K. 

In order to discuss the disordering of the magnetic is- 
land array, we distinguish between the microscopic mag- 
netic moments within the islands (micro-spins, B) and the 
island macro-spins, composed of all n moments within an 
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micro-spins is due to the exchange interaction, with a 
coupling constant J related to the Tc of the magnetic 
material, while the coupling between the islands is mag- 
netostatic. These interactions are long ranged but the 
characteristic energy scale can be defined through the 
interaction between macro-spins, 5, on neighboring is- 
lands. The most relevant interaction is between the near- 
est neighbours, coupling spins along rows or columns of 
the array. We define this coupling as = Kqw? ^ where 



Kq is of order n? and m? 



Introducing the thermally averaged quantity m? im- 
plies that thermal fluctuations of the macro and micro 
degrees of freedom are partially decoupled, allowing the 
two sets to order at different temperatures. This assump- 
tion could break down for smaller units, in which case 
the two ordering processes could be renormalized into a 
single transition, but it should be valid for the islands 
studied here, consisting of ~ 10^ micro-spins. Order- 
ing of the macro-spins due to magnetostatic interactions 
should therefore occur at a temperature Tm < Tc, as 
even in the limit Kq ^ J, the term w? ensures that 
K < J at Tc- Thus we expect disordering of the macro- 
spins at Tm and a range of temperatures Tm < T < Tc, 
in which the islands are decoupled. 

The macro-spin *S is confined to lie along the long axis 
of the island through the magnetic shape anisotropy Ks, 
which minimizes the collective dipolar energy of an is- 
land for the moment in one of the two configurations with 
Ising symmetry. The energy barrier for reversing the di- 
rection of S can be approximated by Er(T) = Koy ^ with 
Ks{T) = fioD{Mom)'^ and is temperature dependent 
through m^. Here jhq is the permeability of free space, 
D is the demagnetization factor of the island, and V 
its volume. The active magnetic material is assumed to 
have a thickness of 1 nm [ ^] when estimating D. The 
macro-spin reversal time scale can be assumed to follow 
a Neel-Arrhenius law r = roe^^/^^^-^^ where /cb is the 
Boltzmann constant and for inter-island magnetic fluctu- 
ations to occur, this barrier must be thermally accessible. 

The temperature dependence of the effective moment 
of the islands is captured by using = (1 — T/Tc)^^, 
with Tc = 230 K, and (3 = 1/3, which matches the tem- 
perature variation of the remanent magnetization of the 
continuous film shown in Fig. 2 to a good enough approx- 
imation for our purposes here. This leads to = E^m? ^ 
with Eo/kB ~ 700 K and Er/{kBT) = 1 at 195 K, which 
ensures accessible thermal fluctuations for macro-spins 
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below Tq. 

The possible vertex configurations for artificial square 
spin ice are summarized in Fig. 1(d). Two dimensional 
square ice does not have the same narrow band of low 
energy states as three dimensional spin ice, as the long 
range part of the dipolar interaction is less well screened 
than in the three dimensional counterpart [ ]. As a 
result, the band of Pauling states cannot be considered 
quasi-degenerate at any level of approximation. The low- 
est energy configurations are made up of type I vertices, 
leading to an ordered ground state, with staggered mag- 
netic moments. The type I vertices, which carry no mo- 
ment, are effectively removed during field cycling and 
are therefore irrelevant in the absence of thermal excita- 
tions. The lowest energy states which both satisfy the 
ice rules and carry a magnetic moment are the subset of 
type II vertices. Hence the simplest model, that captures 
the essential physics observed in this multi-scale problem 
is a sixteen-vertex model in which the type I states are 
raised in energy, leaving a four- vertex (type II) manifold 
of ground states. 

Within this scenario the application of a saturation 
field along the [10] direction would result in a partially 
lifted degeneracy at low temperature, in favor of the en- 
semble of states shown in Fig. 1(b). Applying the field 
along the [11] direction, followed by a return to rema- 
nence would lead to the selection of the fully ordered 
state illustrated in Fig. 1(c). To capture the tempera- 
ture dependence of the excitations we can therefore com- 
pare the ratio of the remanent and saturation magneti- 
zation, Mrem/^sat, for thcsc two dircctious. The order 
parameter of the islands, m^, cancels out, leaving the 
ordering behavior of the macro-spins S only. The result- 
ing ratio should therefore be 1/2 for the [10] direction 
and 1/a/2 for the [11] direction. As seen in Fig. 3), the 
measured ratio falls on the prediction up to 125 K, illus- 
trating collective behavior of the macro spins, with the 
single magnetic domain structure illustrated in Fig. l(b, 
c). Remarkably, the four vertex behavior, in which ver- 
tex defects can be neglected, remains up to temperatures 
where thermal fluctuations have reduced the squared or- 
der parameter of the islands, m^, to less than one half of 
its zero temperature value. In the case of the [11] field, 
the data falls below the 1/a/2 value at 125 K, while that 
for the [10] field stays on the line up to a slightly higher 
temperature. The magnetization, for the field applied in 
both directions, finally collapses as the ordering temper- 
ature of the material, Tq = 230 K, is approached. 

We performed Monte Carlo simulations on a system 
with N = 45000 macro spins, to model the experimental 
results. A sixteen-vertex ice model of N Ising spins S 
was used, arranged on a square lattice. To mimic the 
experimental conditions the energy of type I vertices was 
raised above those of type II and for simplicity we set 
them equal in energy to the type IV vertices, giving a 
single energy scale, em — e// = e/y — em = K. This 
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FIG. 3. The ratio of the remanent to saturation magnetiza- 
tion for the array. At low temperatures the ratios are 0.5 and 
0.707 for the [10] and [11] directions. This is consistent with 
an ordering of type II A + B and type III A + B vertices 
as defined in Fig. 1(d). Above 125 K, the macro-spin order 
diminishes ('pre-melts') which is well below the Curie tem- 
perature of the material (Tc = 230 K). The solid black lines 
are numerical simulations for the 16-vertex model, described 
in the text, that account for the thermally induced dynamics 
of the macro-spins. 
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FIG. 4. Illustration of the thermal collapse of chain-like cor- 
relations between macro-spins. The array is approximated by 
an assembly of finite length Ising chains that order below a 
temperature, Tm- The points are numerical data from Monte 
Carlo simulations and the lines the analytical solution for the 
magnetization of an Ising spin chain of length L — 150. The 
extent of this 'finite size decoupling' regime can be varied by 
changing the coupling strength. Inset: The relation between 
the spin coupling constant and the pre-melting temperature 
Tm, extracted from the Monte Carlo simulations. The finite 
size decoupling regime occurs for any value of the coupling 
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model is equivalent to a grill of interpenetrating Ising 
spin chains coupled uniquely through the first neighbours 
along the chains [see Fig. 1(d)]. The ground state is 
thus an = 2^^ degenerate manifold of type II states 
with ferromagnetic ordering along chains (rows) of length 
L = YiVy2 but with disorder among the chains. The pa- 
rameter K defines the effective coupling of macro-spins 
within a chain and its temperature dependence is ob- 
tained by scaling with m^. As expected, there is no phase 
transition to an ordered state at finite temperature and in 
zero field the system remains disordered to T = 0. How- 
ever the model does have remarkable finite size effects 
related to the zero temperature critical point of the one 
dimensional Ising model [ , ] . This effect is illustrated 
in Fig. 4, which shows the temperature evolution of the 
magnetic moment for macro spins on a single row (or col- 

L 

umn) of spins, Msim = for different values of 

Kq/Tq. In each case the data are compared with Mana, 
calculated for a one dimensional Ising chain of L macro- 
spins [ ], with a coupling constant K = K^m? ^ which is 
exact to leading order in 1/L. The Monte Carlo results 
are excellently reproduced by the analytic expression for 
all values of Kq/Tq^ with small differences coming from 
corrections to scaling [ ] . The chain moment approaches 
saturation below a crossover temperature, which we can 
interpret as Tm when modeling the array. It scales loga- 
rithmically to zero with system size, thus for any reason- 
able system size, there is a finite region of temperatures 
for which each chain of the vertex model is ordered. 

This finite size chain ordering is at the origin of the ob- 
served pre-melting as it ensures that the field hysteresis 
protocol, leading to Mrem, will order the chains giving a 
non-zero total moment for the vertex model below this 
crossover temperature scale. For a field along the [10] 
direction, the chains parallel to the field will be ordered 
but those perpendicular will be free, while for the field 
along [11] all chains should be ordered, as in the experi- 
ment. This field induced ordering is confirmed in Fig. 3, 
where we compare calculated and experimental remanent 
to saturation magnetization ratio values. The numerical 
data was taken from simulations using two fitted param- 
eters: Kq ^ 530 K, corresponding to Kq/Tq = 2.3 (see 
Fig. 4), and a remanent applied field 5o = 1-8 x 10~^ T. 
The simulation and experimental data are in good agree- 
ment up to around 200 K, above which the experimental 
data falls more rapidly towards zero. This could arise 
due to the development of magnetic fluctuations perpen- 
dicular to the Ising spin axis or thermal macro-spin flips 
becoming more favorable. These are not taken into ac- 
count numerically. 

The best flt value of Kq compares favorably with calcu- 
lated values: the direct interaction of the island magnetic 
poles (treated as point charges) yields Kq ^ 530 K. In 
addition, this value of Kq yields a crossover temperature. 



Tm ^ 160iC which is compatible with the experimental 
data. The general agreement between experiment and 
simulation and the closeness of the fltted parameters to 
expected values therefore suggests that the sixteen vertex 
model proposed here, is indeed an accurate description 
of the experimental data. 

In both experiment and simulation, the crossover from 
the ordered spin ice regime to the pre-melting regime be- 
gins at a lower temperature for the [11] fleld direction 
than for the [10] direction. Pre-melting begins as vertex 
excitations disobeying the ice rules occur, breaking up 
the correlated chains of macro-spins. In such an excita- 
tion, a vertex of type II becomes one of type III through 
the flipping of the direction of a macro-spin. In the case 
of remanent magnetization along the [10] direction, there 
are vertex excitations that preserve the magnetic mo- 
ment, for example type 11 B ^ type III ^4 + 5, as 
illustrated in Fig. 1(d). However, for Mrem along [11] all 
excitations are moment reducing. This suggests that the 
magnetization ratio will remain near the upper bound 
in the flrst stages of pre-melting for ordering along the 
[10] direction while the ratio will be reduced from this 
bound as soon as defects appear, for ordering along the 
[11] direction. These defects carry accumulations of effec- 
tive magnetic charge and in the case of three dimensional 
spin ice, have been successfully described using the lan- 
guage of emergent magnetic monopoles [32-34]. In this 
case, the term 'monopole' should be used with great cau- 
tion, as long ranged interactions are less well screened for 
four-fold symmetry in a plane, than in spin ice. Here we 
have demonstrated the possibility to measure the order 
to disorder transition in artiflcial magnetic spin ice. The 
possibility of observing emergent monopoles is therefore 
conceivable, following the general approach which we de- 
scribe in the design of spin ice arrays. 
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